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Abstract. The massless sunrise diagram with an arbitrary number of loops is 
calculated in a simple but formal manner. The result is then verified by rigorous 
mathematical treatment. Pitfalls in the calculation with distributions are highlighted 
fvj ' and explained. The result displays the high energy behaviour of the massive sunrise 

' diagrams, whose calculation is involved already for the two-loop case. 
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1. Introduction 



^ ', The Feynman propagator of a free massless scalar quantum field $(x) 

\^ ' fulfilling the wave equation 

o' 

n$(x) = a^5^$(x) = (1) 

o ■ is given in configuration space by the vacuum expectation value 

n , r fj'il. p-ikx 

AHx) = -.(0|n*W*(0))|0)=/p^^,^ 

V^ i i I i I I ., o, 

U ■ = T^—^ — ^ = T^P^ - ir^ix), 2 

<^ . Att^x^-iQ 47r2 x^ Att ^ ' ^ ' 

where T is the time-ordering operator, P denotes principal value reg- 
ularization and 5 is the one-dimensional Dirac distribution depending 
on x^ = x^x^^ = (x°)^ — (x^)^ — (x^)^ — (x^)^ = Xq — xP'. Formally, one 
could write for the massless (n — l)-loop sunrise diagram 



I (fx[^F{xTe^^'' = Tn{k) 
Ski d*kn-i 1 1 



(27r)4 ■ ■ ■ (27r)4 kl + i^kl + iG" ' {k - ki - k2 - . . . - kn-if + iO ' 

(3) 

<>/ [ © 2008 Springer Science+Business Media, Inc. Manufactured in The Netherlands. 
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Figure 1. The (n-l)-loop sunrise diagram. 



This expression above is both ultraviolet and infrared divergent for 
n > 2, consequently the corresponding expressions in eq. Q are all ill- 
defined. Therefore, [Aj7'(3;)]" has to be regularized in order to obtain a 
mathematically well-defined tempered distribution in the dual Schwartz 
space 5'(]R,^). The regularization can be performed by one of the many 
well-known methods ^ 121 13 ISl El ; which all lead to the same result. 
For n = 2 one obtains 



Tf{k) = n 



d^ki 



1 



1 



log 



A;2-hi0 



(27r)4fc2 + ^0 {k-kiY + iQ 
\og\k'^/fiW 



-9(^2), (4) 



4(27r)2 ^\^ ^\ ) 4(27r)2" °' ""'^' 16^ 

where TZ denotes the regularization procedure with renormalization 
scale ^-ji and Q is the Heaviside distribution. Obviously, the regularized 
expression for the one-loop diagram T^{k) is defined up to a constant 
im momentum space or up to a local distribution in configuration 
space ~ 6^^'{x). The scaling symmetry of the Feynman propagator 
Ap{Xx) = X^^Ap{x), AgH, is spontaneously broken by the regulariza- 
tion procedure, such that a corresponding scaling law r2 (A/c) = T2{k) 
does not hold as one might expect naively from the formal definition 
oiT2{k) ineq. ©. 



2. Formal approach 

In order to obtain the analytic result for an arbitrary multi-loop di- 
agram, a purely formal calculation is presented here first. One easily 
derives formally for n ^ 1 

1 4n(n - 1) 



D 



(x2 - iO)" (x2 - iO)"+l ' 



(5) 
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and therefore 

1^-(^>1" = - 4-^(2,p(-4-l)Kn-2)! °""l^^'-''' • <"> 

Since the d'Alembert operator n goes over into —k'^ in momentum 
space, one obtains 

(_^\n(h.2\n~2 "-2 

^n{k) = ,, i ,;;^ ... -r?(A:) + Vcf(fc2)\ n>2, 

"^ ^ 4"-2(27r)2(""2)(n-l)!(n-2)! ^^ ~ 

(7) 
if the assumption is made that the formal result eq. © is also valid for 
the regularized expressions. The additional terms involving the Fourier 
transform of (covariant derivatives of) the (5-distribution in configu- 
ration space reflect the fact that the regularized Lorentz invariant 
product of Feynman propagators is only defined up to local terms given 
by derivatives of the 5-distribution, whose order can be restricted by 
the requirement that the scaling behaviour of ((TJ corresponds to the 
artificial degree of divergence uj = 2(n — 2) of the (n — l)-loop sunrise 
diagram. Further restrictions follow from physical symmetry require- 
ments, depending on the underlying theory in which the diagram plays 
a role (like in perturbative quantum gravity). As an example, for the 
2-loop diagram (n = 3) one has the general result 



'32(2vr)4^°^l ;.| 



r?W = -^JTTTTTTu l°g ^ +c?- (8) 



The simple result given by eq. ((TJ is indeed correct and similar formal 
calculations as presented above can be found in many textbooks (see, 
e.g., [7]). The massive two-loop case of the sunrise diagram which is 
already involved from the computational point of view has been studied 
in detail in numerous papers (see e.g. |H1 1^ imi ITTl IT^ O ITU [TH | ITKl ITH l 
^^ and references therein), and to a certain extent also for the massive 
multi-loop case fT^ ITHl QUI m] . 

3. Properties of tempered distributions 

A rigorous derivation of ((TJ is presented below. For the forthcoming 
discussion, it is advantageous to review some of the useful properties 
of basic tempered distributions appearing in quantum field theory. A 
scalar neutral field $ can be decomposed into a negative and positive 
frequency part according to 

a(fc)e-^^^ + at(fc)e^^^ . (9) 
•/(27r)32|A?P 
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The commutation relations for such a field are given by the positive 
and negative frequency Jordan-Pauli distributions 

A±(x) = -i[^^{x), $^(0)] = -i(0|[$^(x), $^(0)]|0) , (10) 

which have the Fourier transforms 

A^{k) = f (fx A±(x)e^'^^ = T(27ri) Q{±k^)5{k'^). (11) 

The fact that the commutator 

[$(x), $(0)] = iA+(j;) + iA"(x) =: iA(x) (12) 

vanishes for spacelike arguments x^ < due to the requirement of 
microcausality, leads to the important property that the Jordan-Pauli 
distribution A has causal support, i.e. it vanishes outside the closed 
forward and backward lightcone such that 

suppA(x) CF^UF^, V^ = {x\x^ >Q,±x^ >{)} (13) 

in the sense of distributions. A further crucial observation is the fact 
that one can introduce the retarded propagator A^'^*(x) which coincides 
with A(x) on V —{0}, i.e. (A*"^*, tp) = (A, tp) holds for all test functions 
in the Schwartz space (/7G5(]R,^) with support suppt^ C E,^ — y . One 
might be tempted to write down in configuration space 

/y^\x) = e(x°)A(x), (14) 

and to transform this expression into a convolution in momentum space 

A^^*(A;) = j ^ K{p)Q{k - p) , A{p) = -(2^z)sgn(p0)5(p2) . 

(15) 
The Heaviside distribution @{x^) could be replaced by @{vx) with 
an arbitrary vector v E V~^ = V — dV inside the open forward 
lightcone. The Fourier transform of the Heaviside distribution Q{aP) 
can be calculated easily 

G(A;) = IhnJ d'xe{x'')e-''^°e"'<^^"-"^^ = ^!^6^'\k). (16) 

For the special case where k is in the forward lightcone V~^ , one can go 
to a Lorentz frame where k = {k^,0) such that eq. (|15() becomes 

A-i^i-O n^ i f , A(pO) i f ^^ A(tfcO) 
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Hence, for arbitrary k G V^, A'''^* would be given by the dispersion 
relation 

However, the integral in eq. (|18)) is undefined. One can circumvent this 
problem, e.g., by introducing a mass for the field $ as a regulator for 
the A-distribution. For the massive Jordan-Pauli distribution Am{x) 
one gets due to the (5-distribution in 

Am{k) = -(27ri)sgn(A;°)(5(fc2 - m^) (19) 

the expression (kGV^) 

sgn{tk^)5{t^k'^ -rn^) 



A-(t) = /<i<5e^ 



t + iO 



J^i^-^)-K>^ib)\^i_ ,,.^„.,, (,„, 



2^/Wm{l - t + iO) A;2 



m 



As a special case of the edge of the wedge theorem ^U it is known 
that the Fourier transform of the retarded distribution A'''^*(A;) is the 
boundary value of an analytic function r(z), regular in T"*" := H^+iF"'". 
This way one obtains from r(z) = l/(-z^ — m?), z(^T~^, and m — > 

The analytic expression for the Feynman propagator is then recovered, 
which coincides with 

A^^\k) = AHk) + A-ik) = ^^ (22) 

for k £ V~^ . The singular part of the Feynman propagator is given by 

lm{AF{k)) = -7r6{P) (23) 

and can be deduced directly from the causal Jordan-Pauli distribution 
A in an obvious way. 



4. Analysis of causal properties 

The causal properties of the Jordan-Pauli distribution and the retarded 
propagator can be observed in an analogous way in the case of the distri- 
butions describing a multi-loop sunrise diagram. The normally ordered 



multiloop.tex; 1/02/2008; 23:06; p. 5 



6 ANDREAS ASTE 

product of n free field operators : <I>(x)" : is an operator valued distri- 
bution. Therefore, the same is true for the tensor product : $(x)"' :: 
^{uT ■ |23]- The vacuum expectation value of (0| : ^{x)^ :: $(0)" :|0) 
can be extracted by Wick's theorem and leads to i"n!A"'"(x)"' (the 
combinatorial factor n! will be disregarded in the following), which is 
indeed a well-defined tempered distribution. Calculating the retarded 
part of the causal distribution A„(x) which we define for n > 2 by 

A„(x) := i^{A+ixr - (-ir(A-(xr) = 1(0|[: <^{xr :, : ^(0)" :]|0) , 

(24) 
which has causal support, will finally lead to the desired result eq. Q- 
For this purpose, the Fourier transform of A„(x) has to be calculated. 
This can be done inductively. For n = 2, we first evaluate the Fourier 
transform of A^(x) := A^(x)^ given by the convolution 

" (^ / ^'^ 0(-9°)<^('/')0(^° - k')5{{k - qf). (25) 

The integral vanishes outside the closed backward lightcone due to 
Lorentz invariance and the two G-distributions in eq. H25|). Therefore 
one can go to a Lorentz frame where k = {k^ < 0,0), and using the 
abbreviation E = vy^ = 1^1 ^iid exploiting the 5-distribution in eq. 
(f^ one obtains 

A~^{k'<Q^^) = -^^jll5{{k'f + 2k'E)Q{-E-k') = 



(2vr)^ 

jd\q\\q\^{^+m)Q{-m-k')=-^M-k'm{k' 



^ Z,!-.!!-.!^/"^ i_.i\^/ i_.i , r\\ J- ^/ ^'^^^'"0\2\ 



2(27r)A;0 y '^ "^ ' V 2 '^V^^ ''^ ' ' 4(27r) 

(26) 
and for arbitrary k one gets the intermediate result 

^2{k) = -^Q{-k'mk^). (27) 

According to eq. (|24|) we have A2(x) = A^ (x)^ — A+(x)^, and from eq. 
(|11|) follows A^(— x) = — A~*"(x), therefore we obtain for the Fourier 
transform of A2(x) 

A2{k) = ^sgn{k^)Q{e). (28) 

It is obvious that A2{k) is related to the real part of eq. Q). 
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5. Point splitting method in momentum space 

Instead of introducing a mass term as regulator of the distribution 
A2{k), we present here another strategy, namely a point splitting pro- 
cedure in momentum space, which also leads to the desired result. As 
was shown in |2^, the massless causal distribution A2(fc) can be 
split in analogy to eq. ((TH|) by computing a retarded part A2^gj(fe) in 
momentum space of A2(x)e*'^^ in configuration space via the subtracted 
dispersion relation J23] 



Al retik) = - [ dt ^ ^2(^^ + g) (29) 

2,ret\ ) 2tt J {t - iO){l - t + iO) ^ ' 

for k — q^V~^ , g^ < and performing a limit 

A^/*(fc) = lim[A|,,,(fc) + c,], (30) 

where Cq is a diverging q-dependent constant (i.e., Cq is independent of 
k) which can be chosen in such a way that the limit in eq. (|,'-{()j) exists. 
Note that A2(fc + q) is the Fourier transform of A2(x)e*^^. Therefore, 
one has to compute {k^V^) 

The roots of {tk + qf = t'^k^ + 2t{kq) + q^ = are 

ti,2 = ^i-kq±VN) (32) 

with A^ = {kqf - k'^q^ > {kqf and ti > 0, ^2 < 0. The integral in ^ 
becomes 

{i2 oo v 

= Y^(- log l*2| + log \t2 - 1| - log |ti| + log |ti - 1| - ivr). (33) 

Since t\^2 -^ for q ^ the logarithms log |ti^2 — 1| vanish in this limit 
and the integral reduces to 

" 1^2 (log l*i*2| + ^^) = Y^ (log ^ - ^^) • (34) 
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Adding Cq = j^ log '-^ to eq. H34() leads after analytic continuation 



to the (covariant) splitting solution 



IT 



Arw = ii-5i°s|-^^i. (35) 



in accordance with the well-known result for V2 given by eq. (|H). 

6. Results for higher loops 

We proceed with the Fourier transform of A~(a;)'^ =: A3 (x) 
A3-W = j^JdUQ{-q')5{q')A^{k - q) 

d\ e(-g°)(5(g2)e(g° - k^)e{{k - qf) (36) 



(fq 



4(27r)4 
For k = {k^< 0, 0) we obtain 



A3-(fe) = -^ / ^ e{-E - k')@{{k'f + 2k' E) 



4(27r 



4(2;r)3 J ^'^"'"^"'1 = "35(2^^ ' ^'^^ 

hence 

^^^''^ = "32(iF^'®^"^°^®^^'^- ^^^^ 

Inserting a factor (A; — (7)^ = (k^)'^ + 2k' E in eq. (|37() . corresponding 
to a replacement A^(A;) ^ A~^]^(A;) in eq. (|5^. leads after a short 
inductive calculation to (n > 2) 

Kik) = n TTT^ — ,,, ,, , e(-fc°)9(fc^) (39) 

«v 7 4"-i(n-l)!(n-2)!(27r)2"-3 ^ ^ v 7 v ; 

or 

(l.2\n-2 

^n{k) = 77 TTT^ ,., ,^ , sgn(A;°)9(fc2). (40) 

^^ 4"-i(n- l)!(n-2)!(27r)2"-3 ^ ^ ; v v \ ) 

Keeping in mind that {k'^)^~'^ in Fourier space corresponds to the local 
operator (— n)"^^ in configuration space, a retarded part of A„ can 
be obtained directly from the retarded solution A2^* calculated above 
for A2 in momentum space, since {k'^)"'~'^ ^2^^ (k) is a retarded part of 
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(A;^)"^^A2(A;). The actual calculation is straightforward and leads after 
analytic continuation of A2^^(k) to the desired result 

'4"-i(27r)2"-2(„ _ !)!(„ _ 2)! ^°Sl ^ 






for n > 2. 



j=0 

(41) 



Acknowledgements 

The author wishes to thank Cyrill von Arx for generating the figure. 
This work was supported by the Swiss National Science Foundation. 



References 

1. Pauli, W., Villars, F.: On the invariant regularization in relativistic quantum 
theory. Rev. Mod. Phys. 21, 434 (1949). 

2. 't Hooft, G., Veltman, B.: Regularization and renormalization of gauge fields. 
Nucl. Phys. B44, 189 (1972). 

3. Schwinger, J.: On gauge invariance and vacuum polarization. Phys. Rev. 82, 
664 (1951). 

4. Scharf, G.: Finite Quantum Electrodynamics. 2nd ed., Springer Verlag, New 
York, 1995. 

5. Epstein H., Glaser V.: The role of locality in perturbation theory. Annales 
Poincare Phys. Theor. A19, 211 (1973). 

6. Aste, A.: The two loop master diagram in the causal approach. Annals Phys. 
257, 158 (1997). 

7. Cheng, T.-P., Li, L.-F.; Gauge Theory of Elementary Particle Physics. Oxford 
University Press, New York, 1984. 

8. Kallen, G., Sabry, A.; Fourth order vacuum polarization. Dan. Mat. Fys. Medd. 
29, no. 17, 1. (1955). 

9. Tkachov, F.V.: A theorem on analytical calculability of four loop renormaliza- 
tion group functions. Phys. Lett. BlOO, 65 (1981); Chetyrkin, K.G., Tkachov, 
F.V.: Integration by parts: The algorithm to calculate beta functions in 4 loops. 
Nucl. Phys. B192, 159 (1981). 

10. Bauberger, S., Berends, F.A., Bohm, M., Buza, M.: Analytical and numerical 
methods for massive two loop selfenergy diagrams. Nucl. Phys. B434, 383 
(1995). 

11. Berends, F.A., Davydychev, A.L, Ussyukina, N.I.: Threshold and pseu- 
dothreshold values of the sunset diagram. Phys. Lett. B426, 95 (1998). 

12. Fleischer, J., Kotikov, A.V., Veretin, O.L.: Analytic two loop results for selfen- 
ergy type and vertex type diagrams with one nonzero mass. Nucl. Phys. B547, 
343 (1999). 

13. Post, P., Tausk, J.B.: The sunset diagram in SU(3) chiral perturbation theory. 
Mod. Phys. Lett. All, 2115 (1996). 



multiloop.tex; 1/02/2008; 23:06; p. 9 



10 ANDREAS ASTE 

14. Davydychev, A.I., Smirnov, V.A.: Threshold expansion of the sunset diagram. 
Nud. Phys. B554, 391 (1999). 

15. Caffo, M., Czyz, H., Remiddi, E.: The threshold expansion of the two loop 
sunrise selfmass meister amplitudes. Nucl. Phys. B611, 503 (2001). 

16. Passarino, G.: An approach toward the numerical evaluation of multiloop 
Feynman diagrams. Nucl. Phys. B619, 257 (2001). 

17. Aste, A., Trautmann, D.: Finite calculation of divergent selfenergy diagrams. 
Can. J. Phys. 81, 1433 (2003). 

18. Laporta, S., Remiddi, E.: The analytic value of a 4-loop sunrise graph in a 
particular kinematical configuration. Nucl. Phys. B704, 349 (2005). 

19. Laporta, S.: High-precision calculation of multiloop Feynman integrals by 
difference equations. Int. J. Mod. Phys. A15, 5087 (2000). 

20. Groote, S., Pivovarov, A. A.: Threshold expansion of Feynman diagrams within 
a configuration space technique. Nucl. Phys. B580, 459 (2000). 

21. Groote, S., Korner, J.G., Pivovarov, A. A: On the evaluation of a certain class 
of Feynman diagrams in x-space: Sunrise-type topologies at any loop order. 
|hep-ph/0506286| 

22. Streater, R.F., Wightman, A.S.: PCT, Spin, Statistics and All That. Benjamin- 
Cummings Publishing Company, 1964. 

23. Constantinescu, F.: Distributions and Their Applications in Physics. Pergamon 
Press, 1980. 

24. Diitsch, M., Hurth, T., Krahe, F., Scharf, G.: Causal construction of Yang-Mills 
theories: II. Nuovo Cim. A107, 375 (1994). 

25. Aste, A., G. Scharf, G., Walther, U.: Power counting degree versus singular 
order in the Schwinger model. Nuovo Cim. Alll, 323 (1998). 



multiloop.tex; 1/02/2008; 23:06; p. 10 



